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STUDIES  OF  THE  EFFECT  OF 
DEPTH  OF  FOCUS  ON  SEISMIC  PULSES 


I.  Objectives 

The  principal  objective  of  this  research  is  to 
improve  the  accuracy  of  focal-depth  determination  as 
a  means  of  distinguishing  between  natural  earthquakes 
and  explosions.  The  research  has  consisted  mainly  of 
•valuating  a  method  of  focal -depth  determination  which 
spends  on  recognition  of  the  sequence  of  vents  which 
arise  from  the  reflection  of  the  initial  ?  pulse  at 
•he  surface,  and  evaluating  the  factors  controlling 
the  accuracy  of  this  measurement. 

A  second  objective  is  to  study  the  effect  of  the 
nature  of  a  boundary  on  seismic  pulses  refracted  along 
it  both  in  scale  models  and  in  the  real  earth. 

II.  Technical  Status  and  Major  Accomplishment,' 

A.  Evaluation  of  the  Watson-Merdler  method  of  focal - 

depth  determination 

An  earlier  technical  report  (AFCRL-64-755)  des¬ 
cribed  the  Watson-Merdler  method  of  focal-depth  deter¬ 
mination  in  detail.  A  summary  of  this  is  in  press  in 
a  VESIAC  conference  report.  In  brief,  the  method 
assumes  that  the  seismogram  consists  of  the  sum  of 
three  signals:  P(t),  the  train  of  pulses  eminating 
from  the  focus  and  travelling  to  a  distant  recorder] 

G(t),  a  similar  train  of  waves  reflected  from  the  earth’s 
surface  with  a  reflection  coefficient  arriving  a 
time  At  behind  P(t);  and  n(t),  any  other  signals 
present,  which  are  considered  to  be  noise. 

S(t)  =  P(t)  +  G(t)  +  n(t) 

=  P(t)  +  RRG(t-At)  +  n(t) 
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It  is  assumed  that  one  of  the  effects  of  Gft)  is  to 
complicate  the  seismogram  bycnanging  the  number  and 
size  of  peaks  and  troughs  which  rise  about  the  back¬ 


ground  noise  level.  As  a  means  of  identifying  G(t) 
in  the  seismogram,  S(t)  is  operated  on  by  a  series  of 
inverse  filters  which  find  an  approximation  of  P(t) 
for  assumed  values  of  RR,  At  and  noise  level  (u).  A 
concentration  ratio,  CR,  is  calculated  for  each  P(t) 
found  by  deconvolution  with  the  inverse  filter  using 


the  formula 


(Z( [P(t)|-u)2)for 

all  1p(  t )  1>  u 

T  2(P(t))‘=i 


where  T  is  the  time  during  which  P(t)  exceeds  u. 

It  is  postulated  that  CR  for  a  range  of  At's  close 
to  the  true  time  delay  of  the  surface  reflection,  pP, 
behind  the  direct  compressional  pulse,  P,  from  any 
seismic  event  will  produce  a  CR  which  is  greater  than 
the  concentration  ratio  of  the  original  seismogram, 

CR0* 

Preliminary  tests  showed  that  for  earthquakes  and 
blasts  and  in  scale  models  a  Cj^CRq  was  usually  found 
at  values  of  At  and  RR  close  to  what  would  be  expected. 
It  looked,  therefore,  as  though  the  method  might  provide 
a  simple  means  of  measuring  depth  of  focus,  and, 
therefore,  a  simple  means  of  separating  blasts  from 
all  except  the  shallowest  of  earthquakes. 

Further  tests  showed  that  the  situation  was  not  so 
simple  as  this.  First,  it  was  found  that  there  are 
five  different  causes  which  can  produce  a  value  of 
CR  greater  than  CR0: 

1.  Deconvolution  removes  pP  sequence  from  the  seis¬ 
mogram  (the  solution  sought). 

2.  Deconvolution  removes  sP  sequence  or  some  other 
sequence  than  pP. 
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3.  Deconvolution  enlarges  peaks  on  the  seismogram 
producing  a  larger  maximum  amplitude  in  the 
deconvolved  seismogram  than  in  the  original. 

(C^  is  very  sensitive  to  the  height  of  the 
largest  peak. ) 

4.  Deconvolution  cancels  out  a  later  half  cycle 
or  more  of  a  long  oscillatory  original  seismic 
pulse. 

5.  Deconvolution  phases  random  noise  so  as  to  add 
to  or  subtract  from  the  principal  peaks  of  the 
seismic  pulses. 

As  a  result  of  this,  for  any  seismcfgram  there  are  many 
cases  of  large  CD,  and  the  right  value  must  be  selected 

n 

from  the  set.  This  works  out  to  be  a  process  of  testing 
the  deconvolved  seismogram,  P( t ) ,  to  determine  the  cause 
of  each  large  value  of  C_.  The  third  cause  in  the  above 

n 

list  appears  to  be  the  major  source  of  false  indications. 
Rules  for  simple  identification  are  being  sought. 

The  second  problem  encountered  was  that  in  some 
cases,  no  large  value  of  was  found  where  the  true 
solution  was  expected.  There  are  several  possible 
causes  for  this.  First,  if  is  small,  the  ghost  pulse 
will  be  lost  in  the  noise.  Figure  1  shows  the  expectable 
value  of  Rjj  as  a  function  of  distance  from  the  epicenter 
for  an  assumed  focal  depth  of  10  km  and  the  indicated 
velocity  structure.  Experience  with  the  method  suggests 
that  if  |rh|<0.3,  it  is  unlikely  that  the  ghost  can  be 
regularly  found  in  normal  seismograms.  This  means  that 
only  seismograms  recorded  at  distances  of  22°  and  beyond 
will  be  dependable.  This  limits  the  method  to  events 
large  enough  to  be  recorded  beyond  22°. 

At  the  other  end  of  the  distance  range,  the 
principal  seismic  pulse  recorded  beyond  103°  has  been 
transmitted  thru  the  earth's  core.  Wave  paths  thru  the 
core  are  complicated,  and  for  most  core  distances 
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more  than  one  pulse  is  recorded.  With  more  than  one 
P  pulse  present  there  is  danger  of  finding  At's  due 
to  these  extra  pulses.  For  this  reason,  little  work 
has  been  done  and  none  is  planned  on  pulses  penetrat¬ 
ing  into  the  core. 

Another  problem  is  that  unless  the  length  of  the 
filter  used  in  deconvolving  the  seismogram  is  several 
times  as  long  as  the  seismgram  itself,  the  decon¬ 
volution  is  inaccurate  and  the  concentration  ratios 
found  may  fail  to  indentify  meaningful  solutions. 

This  was  not  realized  until  recently.  As  a  result 
many  of  the  tests  made  for  At^  4  seconds  are  undepend¬ 
able  and  will  have  to  be  redone. 

A  revealing  series  of  tests  of  the  accuracy  of  the 
method  was  made  by  testing  pulses  consisting  of  one 
or  more  cycles  of  a  sine  wave  with  and  without  a  period 
of  quiet  following  the  sine  waves  (Table  I).  It  was 
from  these  tests  that  it  was  first  learned  that  decon¬ 
volution  could  give  large  values  of  CR  where  there  was 
peak-enhancement,  and  that  when  At  was  small,  one  had 
to  be  careful  that  the  deconvolution  filter,  which  is 
proportional  in  length  to  At  in  the  computer  program 
used,  was  long  enough. 

Table  II  summarizes  the  tests  made  to  date  on 
nuclear  blasts.  The  tests  made  for  recordings  at  less 
than  3*000  km  were  early  tests,  often  unsystematic,  and 
the  results  are  of  uncertain  value.  We  have  tested  all 
values  of  At  in  the  ranges  given  in  the  table  in  0,1 
sec  intervals  unless  otherwise  stated  in  the  table.  It 
is  necessary  to  test  at  intervals  which  are  short 
compared  to  the  predominant  frequency  in  the  seismograms 
or  cases  where  CR  is  large  may  be  missed.  The  digitiz¬ 
ing  interval  used  for  the  seismograms  was  0,1  sec  in  all 
cases,  though  data  at  half  this  spacing  is  available. 
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The  quantity,  R,  referred  to  in  Table  I,  is  the 
noise-to-signal  power  ratio  assumed  for  the  seis¬ 
mogram.  This  is  a  factor  on  which  the  inverse  filter 
which  does  the  deconvolution  depends.  The  pulse 
length  tested  is  diortened  for  short  At's  to  keep  the 
inverse  filter  reasonably  short.  This  was  not  done 
in  the  earliest  experiments.  As  a  result,  there  is 
overlap  of  the  At  ranges  tested  as  summarized  in 
Table  II.  It  is  believed  that  most  of  all  of  the 
cases  of  increased  C^0  reported  in  Table  II  will  be 
explainable  by  one  of  the  reasons  listed  above. 

Which  case  corresponds  to  which  cause  is  the  subject 
of  current  study. 

Fourier-integral  frequency  spectra  for  24  of  the 
seismograms  of  Table  II  have  been  calculated.  The 
peak  frequency  and  the  range  for  which  amplitude  is 
within  six  decibels  of  the  peak  are  given  in  Table  III. 
It  is  reasonable  to  expect  that  testing  and  sampling 
intervals  should  be  shorter  than  one-tenth  the  period 
of  the  peak  frequency.  No  thorough  test  of  how  close 
the  sampling  interval  should  be  has  been  made. 

To  evaluate  the  Watson-Merdler  method  it  is 
necessary  to  compare  its  depth  predictions  for  blasts 
with  its  predictions  for  earthquakes.  To  test  a  given 
earthquake  record  thoroughly  for  all  possible  amplitude 
ratios  (RH=pP/P)  and  delay  times  (At)  within  reason¬ 
able  limits  involves’  a  considerable  expenditure  of  time 
and  money.  Initial  detailed  analysis  of  the  Ecuador 
earthquake  of  May  10,  1963*  has  been  restricted  to  a 
reasonably  small  number  of  amplitude  ratios  and  time 
delays  by  considering  the  following  theoretical  model. 
If  one  assumes  that  the  earthquake  can  be  represented 
by  a  point  source  with  energy  radiation  equal  in  all 
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directions,  then  the  theoretical  amplitude  ratio  can 
be  calculated.  This  assumption  is  not  generally  good 
as  energy  radiation  is  known  to  be  dependent  on 
azimuth.  A  knowledge  of  the  initial  radiation  nattern 
(obtainable  from  fault-plane  studies)  would  serve  to 
improve  the  model  but  is  not  available  Huring  this 
stage  of  the  work.  Knowledge  of  the  approximate  depth 
of  the  disturbance  can  be  used  to  restrict  the  range 
of  delay  times  to  be  tested.  Finally,  a  measure  of 
the  noise  level  of  the  seismogram  allows  one  to  pick 
the  set  of  inverse  filters  designed  for  that  partic¬ 
ular  noise-to-signal  power  ratio  and  to  limit  the 
number  of  noise -energy  levels  examined  in  the  actual 
test. 

To  provide  the  necessary  information  for  limiting 
the  length  of  analysis,  programs  were  written  to  give 
the  following: 

(1)  The  angle  of  incidence  of  p  at  the  free  surface 
as  a  function  of  epicentral  distance  for 

10  km  focal  depth  (Figure  2).  The  velocity 
distribution  assumed  for  this  calculation  is 
shown  on  the  Figure . 

(2)  The  amplitude  ratio,  RR,  of  pP/p  as  a  function 
of  the  angle  of  incidence  of  p  at  the  free 
surface  for  various  values  of  Poisson's  ratio. 
The  curve  for  a  Poisson's  ratio  of  0.25  is 
shown  in  Figure  3. 

(3)  Finally,  the  amplitude  ratio  of  pP/p  as  a 
function  of  epicentral  distance  was  detei - 
mined  and  is  shown  in  Figure  1  for  a  focal 
depth  of  ten  kilometers. 

The  pP  wave  is  reversed  in  polarity  compared  with  P 
for  epicentral  distances  greater  than  about  19°,  so 
negative  filters  are  called  for.  It  must  be 
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recognized,  however,  chat  the  focal  mechanism  may  be 
such  that  both  pP  and  P  have  the  same  polarity. 

Preliminary  tests  have  been  carried  out  on  twelve 
records  for  the  10  May  1963  Ecuador  earthquake  using 
the  results  from  the  theoretical  calculations  based 
on  a  point-source  model,  and  some  tests  have  been  made 
on  eight  additional  seismograms.  Pertinent  results 
are  presented  in  Tables  IV  and  V.  The  first  five 
records  listed  in  Table  IV  were  tested  for  delay  times 
between  2.5  and  12.4  sec  corresponding  to  a  depth  of 
focus  in  the  range  of  about  6  to  49  km.  As  no  signif- 
icar  t  answers  were  developed  for  times  greater  than 
6  sec.  Initial  tests  on  the  remaining  seven  records 
were  only  carried  out  to  8  sec  (about  30  km).  Depths 
less  than  6  km  have  not  been  tested,  as  the  present 
set  of  inverse  filters  are  not  designed  properly  for 
che  small  delay  times.  Table  IV  shows  the  theoretical 
amplitude  ratio  and  measured  noise -to-signal  power 
ratio  for  each  record  and  the  filters  used  in  the 
tests.  In  all  except  two  cases  there  was  more  than  one 
case  where  CR  exceeded  CRQ.  Only  the  solutions  for 
the  largest  C^'s  are  listed  in  Table  IV.  All  solutions 
are  shown  in  Figure  4. 

No  consistent  depth  persists  throughout  Table  IV. 
Values  range  from  8  to  20  km.  The  average  value  is 
12  km  (compared  to  30  km  reported  by  the  Coast  and 
Geodetic  Survey) .  In  Figure  4,  on  the  other  hand, 
there  is  a  solution  between  8  and  12  km  in  every  case 
except  where  no  solution  was  found  at  any  depth. 

Additional  test  results  using  filters  with  a  high 
noise-to-signal  power  ratio  of  0.25  are  included  in 
Table  V.  The  amplitude  ratios  giving  significant 
answers  are  in  general  higher  than  those  listed  in 
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Table  IV.  The  spread  in  answers  is  smaller  and 
there  appears  to  be  a  preferred  depth  of  about 
.12  km.  Causes  of  the  numerous  answers  for  eacn 
record  are  being  sought  as  we1!  as  a  means  for 
picking  the  "best”  answer. 

In  addition  to  those  of  the  10  May  1963 
earthquake,  nineteen  other  seismograms  have  been 
processed  thru  the  Watson-Merdler  program.  These 
are  listed  in  Table  VI.  The  thoroughness  of 
these  tests  varies  from  scattered  trial  solutions 
to  systematic  tests  following  the  same  systems 
used  for  the  nuclear  blasts.  These  tests  have 
been  run  for  a  variety  of  purposes  such  as  for 
comparison  with  nuclear  blasts  and  as  checks  to 
see  if  any  peculiar  solutions  would  result  from 
various  features  of  the  records.  They  have  led 
to  no  conclusions  not  covered  by  the  previous 
discussions  of  nuclear  blasts  and  the  10  May  1963 
Ecuador  earthquake. 

B.  Factors  controlling  the  frequency  spectra  of 

seismograms 

In  order  to  get  a  better  understanding  of  tne 
factors  which  control  the  frequency  spectrum  of  a 
seismogram,  a  study  has  been  made  of  the  earthquakes 
of  13  April  1963  in  Peru  and  10  May  1963  in  Ecuador. 
As  many  seismograms  as  possible  of  the  beginning  of 
this  earthquake  were  obtained  from  the  U.  S.  Coast 
and  Geodetic  Survey  and  the  Dominion  Observatory  in 
Ottawa,  Canada.  As  most  of  the  effort  was  directed 
at  the  Peru  earthquake  it  will  be  the  only  one  of 
the  two  for  which  detailed  results  will  be  summarized 
here . 

Of  the  paper  seismograms  obtained  for  the  Peru 
earthquake,  only  ten  (five  from  USCGS  and  five  from 
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Dominion  Observatory)  were  sufficiently,  clearly 
and  continuously  recorded  that  they  could  be  used 
( see  Howell,  1965  for  discussion  of  problems  of 
digitizing  paper  seismogram).  In  each  case  the 
seismogram  was  digitized  at  tenth-second  intervals 
for  the  first  10.0  ±  0.4  sec  of  the  seismogram.  Twenty- 
seven  additional  digitized  seismograms  were  obtained 
from  the  United  Electrodynamics  data  center  in 
Alexandria,  Virginia.  Fourier-integral  TVequency 
spectra  were  calculated  and  a  correction  for  station 
sensitivity  applied  over  the  frequency  range  0.4  to 
2.5  cps  for  each  of  these  37  seismograms.  The  range 
of  epicentral  distances  was  from  4411  km  to  11,326 
km.  Five  other  Vela  stations  located  in  the  shadow 
zone  or  further  away  also  gave  excellent  recordings 
but  were  not  included  in  this  study  of  the  P-puise. 

Samples  of  the  corrected  frequency  spectra  are 
shown  in  Figure  5.  Except  for  a  general  decree  ie  in 
energy  with  increase  in  frequency,  the  spectra  exhibit 
no  systematic  pattern.  Maxima  and  minima  occur  in  all 
cases,  but  their  positions  and  heights  are  different 
in  nearly  every  case. 

A  report  by  W.  J.  Hannon  ( 1964a  and  B)  led  to  a 
possible  explanation  of  this  extreme  variability. 

Hannon  developed  a  program  for  computing  the  effect 
of  a  layered  crust  on  the  interference  pattern  of 
continuous  plane  waves  being  passed  thru  it  and 
reflected  at  the  surface.  The  spectrum  observed  at 
the  surface-is  the  result  of  interference  by  the  whole 
family  of  internally  reflected  multiples  resulting  from 
the  layering.  The  effect  of  these  layers  on  the 
spectrum  is  called  the  transfer  function  of  the  crust. 
The  transfer  function  itself  is  a  sequence  of  maxima 
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and  minima  (Figure  6). 

Hannon's  theory  is  developed  for  continuous 
waves.  A  real  earthquake  is  a  transient  pulse.  If 
only  the  first  ten  seconds  of  a  seismic  arrival  is 
considered,  many  of  the  internal  reflections  contrib¬ 
uting  to  Hannon1 s  theoretical  transfer  function  will 
generally  not  have  had  time  to  arrive.  It  was,  there¬ 
fore,  necessary  to  develop  a  program  for  calculating 
the  transfer  function  for  the  case  of  a  finite  length 
of  pulse 

Truncation  in  the  time  domain  corresponds  to 
multiplication  of  the  signal  by  a  square  window  of 
unit  amplitude  during  the  time  duration  of  the 
considered  signal,  and  equal  to  zero  outside  this 
interval.  The  transform  of  such  a  square  window  is 

^  sine  function  whose  zero-axis  crossing  will  be  at 
ip  cps,  T  being  the  duration  in  seconds. 

Cc  evolution  in  the  frequency  domain  is  equivalent 
to  multiplication  in  the  time  domain.  Let  G(f)  be 
the  transform  for  infinite  time  with  9(f)  its  related 
phase  spectrum.  The  truncated  spectrum 

G(f)  =  G(f)*sinc(f) 

|aT(f)|  =((ReG(f)*sinc(f))2  + 
(ImG(f)*sino(f))2)1/2  (4) 

Proper  caution  must  be  taken  during  the  convolution 
process  to  reduce  as  much  as  possible  the  error  intro¬ 
duced  b>  the  truncation  of  the  sine  function.  This  is 
done  by  extending  the  real  and  imaginery  parts  along 
tl  e  negative  frequency  axis  before  performing  the 
convolution. 

Computer  programs  were  written  to  generate  transfer 
functions  corresponding  to  ten  seconds  in  the  time 
domain  for  angles  of  incidence  at  every  1°  within  the 
range  of  20°  to  40°.  This  range  includes  all  expected 
angles  of  incidence  for  epicentral  distances  involved 
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In  the  study  of  the  37  spectra. 

More  than  20  crustal -structure  models  related 
to  seven  stations  were  studied.  It  was  found  that 
certain  crustal  models  have  a  smooth  transfer 
function  for  the  whole  range  of  angles  of  incidence, 

■  hile  other  models,  over  the  same  range  of  angles, 
present  great  differences  in  the  number,  the  position 
and  the  depth  of  the  minima.  Since  the  crustal  effect 
acts  as  a  filter  on  the  incoming  initial  longitudinal 
wave,  it  modifies  the  frequency  spectrum  of  the  pulse 
greatly. 

The  variation  of  the  truncates  transfer  function 
as  a  function  of  angle  cf  incidence  and  of  layer 
thicknesses  was  carefully  studied  by  testing  nine 
closely  related  models  selected  on  the  basis  of  pub¬ 
lished  crustal-structure  predictions.  It  was  found 
that  the  transfer  function  is  very  sensitive  to  change 
of  angle  of  incidence  and  layer  thickness.  The 
transfer  functions  can  sometimes  introduce  minima  of 
20-30  db,  very  often  of  10  db,  thus  offering  an  adequate 
explanation  for  similar  minima  found  in  earthquake 
spectra. 

An  attempt  was  made  to  find  crustal  models  which 
would  produce  variations  similar  to  those  observed 
in  the  spectra  at  four  stations.  In  two  of  these 
cases,  namely  Mountain  Pine,  Ark.  and  Alert,  Canada, 
a  structure  was  found  whose  transfer  function,  sub¬ 
tracted  from  the  observed  spectrum,  produced  a 
demonstrably  smoother  spectrum  than  the  original. 

It  is  postulated  that  this  model  closely  resembles 
the  real  crust.  Failure  to  find  a  model  which  produced 
an  improvement  for  the  other  two  cases  is  presumably 
due  to  not  having  hit  upon  a  model  which  closely 
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resembles  the  crust  among  the  many  models  tested. 

Finally,  the  averse  spectrum  was  obtained  for 
the  37  stations  (Figure  9)*  It  shows  a  remarkable 
smoothness  in  comparison  with  individual  spectra. 

This  is  interpreted  as  an  indication  that  none  of 
the  minima  is  introduced  by  the  source  function. 
Otherwise  it  would  persist  and  show  up  on  the 
average  spectrum. 

It  is  concluded  that  the  crustal  transfer 
function  is  th'  dominant  cause  of  variations  in 
the  spectrum  from  one  observing  station  to  another. 

C.  Scale -model  studies  of  refractions  along  the 

Mohorovlclc  discontinuity.  * 

Two  scale  models  were  designed  and  fabricated 
as  shown  in  Figure  10.  The  upper  model  (A)  corres¬ 
ponds  to  a  transition  layer  underlying  a  uniform 
crust.  The  lower  model  (B)  corresponds  to  an 
irregular  boundary.  Model  A  has  been  set  up  in  the 
model  laboratory  and  is  under  test. 

In  Model  A  an  attempt  is  being  made  to  measure 
the  transfer  characteristics  of  the  two  high-speed 
interfaces  comprising  an  intermediate  velocity  layer, 
with  the  distance  between  these  interfaces  being  of 
the  same  order  of  magnitude  as  the  px’edominent  wave 
length.  The  frequency  band  of  the  study  is  from  20 
kc  to  200  kc,  which  is  approximately  equivalent  to 
20  cm  to  2  cm  wave  length  in  the  model.  The  irregular 
boundary  at  the  bottom  .f  the  model  is  an  attempt  to 
reduce  reflections  occuring  at  this  surface,  prevent 
the  arrival  of  surface  waves  that  may  travel  around 
the  model  by  this  path,  and  reduce  the  time  that 
energy  is  contained  in  the  model  between  pulses. 
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An  electronics  system  capable  of  providing  a 
short- duration  input  to  the  model  was  assembled. 

The  system  is  able  to  put  out  square  pulses  of  as 
short  a  duration  as  one  microsecond  and  of  150 
volts  amplitude.  Different  preamplifiers  and  even 
the  cascading  of  preamplifiers  was  tried  to  find 
out  how  much  gain  would  be  available  when  low-energy 
refractive  arrivals  are  to  be  measured.  A  high 
signal-to-noise  ratio  is  necessary,  and  precautions 
had  to  be  taken  to  :  nsure  low  pickup  and  interfer¬ 
ence.  This  will  become  critical  at  the  lowest 
energy  levels  anticipated. 

The  velocities  of  compressional  waves  in  the 
three  layers  were  measured  (Table  VII).  In  addition, 
a  calculation  of  the  velocity,  C  ,  in  the  intermediate 
velocity  layer  (copper-steel  lamination)  was  made 
using 

P-.Cni^n  +  Po^do^o 

(5) 


2  _  P1CP1T1  +  P2CP2T2 
CP 


PlTl  +  p2t2 

where  p^,Cp-^,T,  are  density,  compressional  wave 
velocity,  and  th'.ckness  of  material  number  1  and  where 
p2,CP2,T2  are  density*  compressional  wave  velocity, 
and  thickness  of  material  number  2. 

The  shear-wave  velocities  were  calculated  several 
ways,  and  also  roughly  measured,  with  an  approximate 
value  being  needed  to  identify  fully  the  components 
present  in  the  preliminary  measurements. 

Measurements  of  the  absorptive  and  dispersive 
characteristics  of  the  materials  were  made.  This 
data  is  currently  under  study.  Having  measurements 
iur  the  properties  of  the  steel  and  the  laminated 
material,  the  copper  material  characteristics  can  be 


Inferred.  Preliminary  evaluation  of  the  data  suggests 
that  the  three  materials  will  be  low  absorptive  and  non- 
dispersive . 

D.  Attentuation  of  seismic  waves  under  Lake  Superior 
Fourier-integral  frequency  spectra  were  calculated 

for  the  range  4-26.8  cps  for  first  arrivals  in  two 
distance  ranges:  35-82  km  and  158-252  km  for  seismic 
pulses  refracted  under  Lake  Superior.  The  first  case 
is  interpreted  as  a  direct  pulse,  the  second  as  a 
refracted  (head)  wave  in  an  intermediate  (7*00  km/sec) 
layer.  The  absorption  coefficient  in  the  first  case 
was  found  not  to  be  significantly  different  from  zero 
(a  =(-0.010151  -  .00059F)/km,  where  F  is  frequency). 

In  the  second  case,  absorption  was  larger  but  not  signif¬ 
icantly  dependent  on  frequency  (a  =  (0.006620  -  0. 00011F)/km) . 
The  refracted  pulse  was  weaker  in  energy  above  8.8  cps 
than  the  direct  pulse.  This  is  similar  to  the  filtering 
effect  observed  in  the  Maine  seismic  experiment  and  is 
explainable,  using  Nakamura's  theory,  by  a  transition 
layer  0.l6  km  thick.  (Nakamura,  1964;  Nakamura  and 
Howell,  1964).  A  report  summarizing  this  is  in  press 
(Howell,  I966A). 

E.  Publications  and  oral  reports 

In  addition  to  the  reports  required  by  this  contract, 
during  1965  the  following  papers  were  presented  orally 
or  in  writing  thru  the  indicated  media,  thereby  dissemi¬ 
nating  knowledge  gained  thru  this  research. 

1.  Robert  J.  Watson  presented  a  paper  describing  the 
Watson-Merdler  method  of  focal-depth  analysis  before  a 
VESIAC  conference  in  La  Jolla  on  March  23.  The  proceed¬ 
ings  of  this  conference  have  not  yet  been  published,  but 
are  to  be  issued  in  1966. 
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2.  B.  P.  Howell,  Jr.  presented  a  paper  before  the 
Seismological  Society  of  America  discussing  some  of 
the  problems  of  obtaining  digiti2able  seismograms 
using  paper  recordings.  This  paper  has  been  published 
in  Earthquake  Notes  (Howell,  1965). 

3.  B.  P.  Howell,  Jr.  presented  a  paper  before  the 
Seismological  Society  of  America  describing  the 
digitizing  apparatus  developed  in  the  research  and 
used  for  digitizing  seismograms.  This  paper  (Howell, 
1966b)  will  be  published  shortly  and  distributed  as  a 
scientific  report. 

4.  B.  F.  Howell,  Jr.  has  prepared  a  report  of  his 
analysis  of  the  frequency  spectra  of  seismic  refractions 
under  Lake  Superior.  This  is  being  submitted  for  publi¬ 
cation  in  the  Merle  A.  Tuve  Testimonial  Volume  of  the 
Carnegie  Institution  of  Washington. 

III.  PERSONNEL 

The  personnel  engaged  in  this  research  during  the 
past  year  consisted  of  Professors  B.  P.  Howell,  Jr., 

P.  M.  Lavln,  and  R.  J.  Watson  and  six  graduate  assistants: 
Y.  Y.  Cheng,  G.  Leblanc,  J.  L.  Lin,  T.  W.  Novotry,  D.  E. 
Siskind,  and  S.  K.  Yiu.  Professor  Howell  has  been  guiding 
the  studies  of  actual  earthquakes  by  Leblanc,  and  Yiu  and 
the  scale-model  work  of  Siskind.  Professor  Lavin  has 
guided  the  scale  model  experiments  of  Novotny  and  an 
evaluation  of  uhe  Watson-Merdler  method  as  applied  to  the 
10  May  1963  earthquake  by  Cheng.  Professor  Watscn 
directed  Mr.  Lin  in  evaluating  the  Watson-Merdler  method 
of  depth-of-foous  calculation  for  nuclear  blasts  until 
he  resigned  in  July;  Professor  Howell  has  directed  Lin 
since  then.  Professor  Howell  has  been  studying  the 
spectra  of  pulses  refracted  thru  the  earth’s  crust  under 
Lake  Superior. 
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IV.  Problems  Encountered 

We  have  been  plagued  throughout  the  work  by 
slow  delivery  of  parts  and  equipment.  During  part 
of  the  past  month  the  University's  computer  has  been 
closed.  As  a  result,  we  are  behind  schedule  and  have 
had  to  obtain  an  extension  of  time  to  complete  the 
proposed  work. 

V.  Schedule  of  work  and  future  plans 

The  research  planned  for  the  remainder  of  this 
project  consists  primarily  of  three  experiments: 

1.  study  of  the  cause  of  large  values  of  for 
deconvolved  blast  seismograms;  2.  study  of  the 
10  May  1963  earthquake;  and  3.  model  studies  of 
refraction  along  the  Mohorovicic  discontinuity. 

Fourteen  nuclear-blast  seismograms  recorded  at 
large  distances  from  the  epicenter  are  on  hand.  All 
are  being  systematically  tested  using  a  standard 
procedure.  In  all  cases  where  CR  exceeds  CRQ,  the 
deconvolved  seismogram  will  be  compared  with  the 
original  seismogram  to  see  exactly  why  the  Watson- 
Merdler  method  picked  this  case  as  a  potential 
solution.  If  all  potential  solutions  have  obvious 
explanations,  the  Watson-Merdler  method  can  be  used 
as  a  criterion  to  select  blast  records  from  a  group 
of  seismograms  of  unknown  source  nature. 

The  work  on  earthquake  seismograms  wil]  follow 
similar  lines.  An  explanation  will  be  sought  for 
all  large  values  of  for  the  10  May  1963  Ecuador 
earthquake  and  for  a  few  other  individual  shocks. 

For  the  Watson-Merdler  method  to  be  successful,  it 
will  3  necessary  that  natural  earthquakes  commonly 
produce  large  values  which  cannot  be  easily 
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explained  as  being  other  than  the  result  of 
combinations  of  pulses. 

Data  from  the  scale  models  simulating  refractions 
along  the  Mohorovicic  discontinuity  will  be  taken 
during  the  coming  quarter.  Analysis  of  this  data 
should  be  complete  by  June. 

Several  reports  are  under  preparation.  A 
Complete  technical  report  and  one  or  more  scientific 
reports  on  the  cause  of  variations  in  the  frequency 
spectra  of  the  13  April  1963  Peru  earthquake,  including 
the  program  for  deriving  the  truncated  transfer  function 
for  a  layered  crust,  is  being  written.  Also,  a 
scientific  report  on  the  spectra  of  refractions  under 
Lake  Superior  is  nearing  completion.  Some  additional 
research  on  this  latter  problem  continues,  especially 
as  it  relates  to  absorption  of  seismic  waves.  This 
may  lead  to  an  additional  report  later. 

VI.  Actions  desired  of  ARPA 

A  proposal  to  extend  this  research  thru  December, 
1966,  either  with  or  without  additional  funds,  was 
submitted  a  year  ago  and  has  never  been  acted  on 
formally.  It  now  looks  as  though  funds  already  approved 
will  be  adequate  to  complete  what  needs  to  be  done  on 
this  research.  Although  actual  research  should  be 
completed  by  June  of  1966,  it  is  still  desired  the 
contract-terminated  date  should  be  postponed  until 
31  December  1966,  This  will  permit  the  final  reports 
to  be  prepared  more  carefully  and  thoroughly  than  will 
be  possible  if  they  must  be  in  by  June.  Also,  it  will 
permit  preparation  and  submission  of  a  series  of 
scientific  reports  for  publication  in  journals  where 
the  information  they  contain  will  reach  the  largest 
number  of  scientists  and  will  become  a  matter  of  con¬ 


venient  permanent  record. 
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Figure  1  -  Epieentra.l  Distance  (Degrees) 
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Figure  2,  Angle  of  incidence  of  P  vs.  epicentral  distance 
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Figure  4.  Chart  of  cases  where  CR  exceeds  CR0  for  seismograms  listed 
in  Table  4. 


Frequency  in  cps 


Figure  6,  Truncated  transfer  function  for  Ontario,  Angle 
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Figure  10.  The  two  models  being  tested 


